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ABSTRACT 

We investigate the extraction of the rotational energy of a black hole under different accreting 
environment. When the accretion rate is moderate, the accretion disk consists of an outer thin disk 
and an inner advection-dominated accretion flow. In such a combined disk, the outer thin disk can 
sustain a magnetic field with moderate strength at the event horizon, leading to the formation of 
relativistic jets with moderate luminosity and speed via the magnetohrodynamic Penrose process. 
When the accretion rate increases enough, on the other hand, the disk becomes geometrically thin 
near the horizon. In this slim disk, the denser plasmas can sustain a stronger magnetic field than 
that in a combined disk, leading to the formation of jets with greater luminosity and speed via the 
Blandford-Znajek processs. It is discussed that the former jets are associated with the Fanar off- Riley 
(FR) I galaxies and the latter with FR II galaxies. 

Subject headings: accretion, accretion disks — black hole physics — Galaxies: active — magnetic 
fields — MHD 



1. INTRODUCTION 

Remarkable advances have been made in explaining 
the observed spectrum of black hole X-ray binaries 
(BHXBs) and active galactic nuclei (AGNs) with the 
help of theo retical studies on black hole (BH) accre- 
tion disks (jAbramowicz et al.l fl999l : iKato et all 12008( 1. 
Current stable BH accretion disk solutions include the 
advection-dominated accret i on flow, or ADAF for short 
dNaravan fe Yil 11994 119951: INaravan et al.l 11997ft. the 
thin disk (jShakura fc Sun vaev 1973; Novikov & T hornd 
j!973h. and the slim disk ( Abramowicz et al.lll988l 120101 : 
|dowski 2009; S adowski et al.ll2011f ). The accretion disk 
type is a function of the the dimensionless accretion rate, 
m = M/MEdd = L/L-Edd, where M is the accretion rate, 
^Edd = ^Edd/£c 2 is the Eddingtion accretion rate, L^dd 
is the Eddingtion luminosity, £ is the efficiency (£ = 1 
is adopted for later computation) and c is the speed of 
light. In Table 1 we summarize the properties of the 
above disk solutions. 

Th e accretion disk solution is also a function of ra- 
dius (jAbramowicz et aT| [l995: 



]hen et al.lll995f ). A com- 
bined disk which consists an outer thin disk and an inner 
ADAF has been also p roposed to explain observations 
of BHXBs and AGNs (INaravan et al 1 119961 : lEsin et al.l 
[19971 : iFender et al.l 120041 : iTrump et al.ll20Tl . Such disk 
is expected to be form before the entire disk transits from 
an ADAF to a thin disk, because the hot ADAF would 
cools down from outside. In other words, the transition 
radius from outer thin disk to inner ADAF solution de- 
creases as the accretio n rate increases until the entire 
disk become thin disk dHp nma 199i| IKato fc Nakamural 
1998; iManmoto et al. 2000|). For a comprehensive con- 
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sideration of how accretion disk type vary with m, we 
also include the the combined disk in Table 1. 

As a pure relativistic effect, the specific angular mo- 
mentum of a Keplerian rotating object near the BH, £^ 
has a local minimum at the radius of the innermost stable 
circular orbit, i^isco- The flow geometry near the BH 
horizon is determined by how the angular momentum 
distribution of the accret ion disk, £(R) , deviates from 
the Keplerian values (jAbramowicz 1119981 ). If £(R) > £ ms , 
where £ ms is the Keplerian augular momentum at i^isco, 
£(R) must equal the Keplerian value twice (assuming 
£{R) = £k at Ri and i?2, and R2 > Risco > Ri)- In this 
case, a pressure maximum is formed inside the disk (at 
R2) and plasmas are "pushed" towards the BH, forming a 
'pressure-driven' flow and resulting in a 'disk-like' struc- 
ture on the equatorial plane near the horizon. On the 
contrary, if £{R) < ^ ms , the sub-Keplerian plasmas fall 
into the BH because of the nature of their sub-Keplerian 
augular momentum ('viscous-driven'). The plasmas are 
able to fall closer to the axis (i.e., higher latitudes of 
the horizon), forming a 'quasi-spherical' geometry near 
the horizon. These two kinds of accretion flow geometry 
near the horiz on, disk-like and quas i -sphe rical was first 
recognized by lAbramowicz fc Zurekl (jl98l[ ) . 

Figure 1 schematically depicts the geometric and dy- 
namic properties of different accretion disks and the ex- 
traction of the BH rotational energy when different ac- 
cretion flow geometries near the horizon is realized. If 
the magnetosphere near the horizon is filled with plas- 
mas, those plasma can load onto the large-scale field lines 
and r educe the electromag netic extraction of the BH en- 
ergy (Takaha shi et al.lll99Q[ ). In this letter, by assuming 
that a relativistic jet are launched when the BH energy 
is extracted outward, we investigate the formation of BH 
relativistic jets when the BH is surrounded by different 
ty pe of disk s listed in Table 1. This work is an extension 
of lPu et all (|2012f ). hereafter PHC, in which the forma- 
tion of relativistic jets at low accretion rate (e.g. when 
m < 0.1 ) is considered. It is found that, ergospheric jets 
can preferentially take place not only when the BH is sur- 
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TABLE 1 

Summary of Black Hole Disk Properties 



Disk Type 



ADAF 



Combined Disk 



Thin Disk 



Slim Disk 



Accretion Rate 

Cooling 
Process 

Optical Depth 
Disk Geometry 

Viscosity 6 

Infalling 

Process 

Flow 
Geometry near 
the Horizon^ 



rh < 0.01 

advection 
thin 

thick (H ~ R) a 

a ~ 0.1 

viscous-driven (l- iri < f ms ) c 
quasi-spherical 



rh < 0.01 



outer thin disk 
+ inner ADAF 



0.01 < rh < 0.3 

radiation 
thick 

thin (H < R) 

a ~ 0.01 

viscous-driven 

quasi-spherical 



0.3 < rh 

advection 
thick 

thick (H < R) 

a < 1 

pressure-driven (£[ n > f ms ) 
disk-like 



a H is the disk height and is the radial distance to the central BH. 
b a is the parameter in (^-prescription (Shakura & Sunvaev 1973). 

c £i n is the angular momentum of the flow when it falls onto the BH. See also Figure 1. 

d Note that the geometry of the accretion flow near the horizon is determined by the 'infalling process' (see text), therefore it can be 
different from the 'disk geometry', which describes the geometry of an accretion disk far from the horizon. See also Figure 1. 



rounded by a combined disk, as suggested in PHC, but 
also when the BH is surrounded by a slim disk. We refer 
the two types of jets as 'EJC' and 'EJS', respectively (see 
section 3). We conclude that the characteristic accretion 
rate of a EJC and a EJS are rh < 0.01 and rh > 0.01, re- 
spectively, and that the EJS is in general more powerful 
and faster than the EJC for the same central BH mass 
and BH spin. The difference between these two types of 
jets may result in the observed dichotomy of of FR I and 
FR II galaxies. 

2. IDEAL MHD FLOW 

In an axis-symmetry and stationary BH magneto- 
sphere, the total energy of the ideal MHD flow, E, is 
conserved along a specific field line (see, e.g., PHC for 
details). We can separate E into the plasma part (the 
first term) and electromanetic part (the second term), 

where \i is the relativistic specific enthalpy, u t is con- 
variant time component of the four — velocity, dp is the 
angular velocity of the field, r] oc nu r /B r is the particle 
flux per unit flux tube, n is the proper number density, 
u r is the radial component of the four-velocity and 
and B r is the toroidal and radial field observed by dis- 
tant observer. Similarly, the total outward energy flux 
S r has the plasma part, £g m , and the electronmagnetic 
part, £g m , 

£ r = nEu r =£: m + £; iasiaa . (2) 

A plasma is accelerated along the large-scale magnetic 
field lines inwards (or outwards) if the gravitational 
force is larger (or less) than the magnetocentrifugal 
forces. The conservation of £ r near the separation 
region, where the magnetocentrifugal forces balance 



with the gravitational force, requires the solution of 
a powerful outflow (i.e. relativistic jet) to that of an 
inflow with a positive energy flux (S r > 0). As discussed 
in §3.6 of PHC, for the inflow, £ e r m > and ££ lasma < 
are generally expected; that is, the magnetic field lines 
contribute to the extraction of the rotational energy 
of a BH, while the plasma plays an opposite role and 
reduces the extraction. Therefore, an MHD flow can 
extract the BH energy (accordingly a ergospheric jet is 
launched) , only when the magnetosphere is magnetically 
dominated, |££ m | > |££i asma |- The extraction of the 
rotational energy by the MH D flow is described b y 
the 'MHD Penrose process' (|Takahashi et all 119901 ). 
When £p lasma is negligibly small compared to ££ m (i.e., 
a nearly vacuum environment), the force- free limit 
becomes a good approximation. The electromagnetic 
extraction of the rotational energy of a BH (by large- 
scale, hole-threading field lines) is described by th e 
'Blandford-Znajek process' (jBlandford fc Znajeklll977l ). 



3. TWO TYPES OF RELATIVISTIC JET 

In general, the strength of a large-scale magnetic field 
strengtlfl can be parameterized by the ratio, e 2 , of the 
gravitational binding energy of the disk at radius R to 
the large-scale magnetic field energy inside R (see also 

4 We should stress here that the estimation of large-scale field 
should be different from the estimation of the local field, which can 
be computed by the plasma beta and the gas pre ssure inside the 
disk. Numerical simulations, e.x. McKinnev et al. (2012) and the 
references therein, show that the accumulated field strength near 
the BH can even reach to an equipartition level and supports the 
accretion disk. 
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quasi-spherical 
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Fig. 1. — (a)Illustration plot (left) and representative angular mo- 
mentum profile (right) of a ADAF. Dominating force are shown: 
F g , force due to the gravity of the BH; F c , force due to the rotation 
of the disk; F p , force due to the pressure gradient inside the disk. 
The dashed line represents the location of i^isco- The specific an- 
gular momentum of the disk is shown in red solid curve while the 
Keplerian value is shown in dashed bule curve. (b)Same as (a), but 
of a thin disk, (c) Same as (a), but of a slim disk. (d)Extraction 
of BH rotational energy by a large-scale magnetic field under dif- 
fe rent accretion flow geometries near the BH horizon (cf. Figure 2 
of Abramowicz & Zurek (1981)). The field lines are indicated by 
solid lines and the accreting plasma are indicated by the shaded 
region (not to scale). (left)Quasi-spherical flow geometry. Such 
geometry is resulting from a viscous-driven acrretion process, such 
as an ADAF or a thin disk. The electromagnetic extraction of 
the BH rotational energy by the field lines is reduced, because the 
inflowing plasmas along the large-scale field lines bring the rest- 
mass energies onto the BH. For the thin disk case, the depicted 
region here is corresponding to a much inner region than i^isco- 
(right)Disk-like flow geometry. Such geometry is resulting from a 
pressure-driven accretion process, such as a slim disk. 



PHC). That is, 



B(s, R, E(i?)) = £v / 27rGM BH S/ J R 2 



(3) 



where G is the gravitational constant and Mbh the mass 
of the central BH. When e = 1, B reaches a equipartition 
value at the radius R. 

By examining whether the BH rotational energy can 
be extracted by the MHD flow via the MHD Penrose 
process when the BH is surrounded by an ADAF, a com- 
bined disk or a thin disk, PHC concluded that the jet 
is most likely launched when the surrounded disk is of 
a combined disk type among all the cases. Because the 
surface density of the outer thin disk is much larger than 
that of the inner ADAF for a combined disk, the inner 
ADAF can therefore 'receive' more large-scale fields from 
the outer thin disk than what a single ADAF would cap- 
ture (see equation [3]). Therefore, the quasi-spherical 
plasma provided by the inner ADAF of a combined disk 
can be loaded on relatively stronger field lines. Relativis- 
tic jets can be finally initiated when the magnetosphere 
becomes magnetically dominated. This situation can be 
realized when the transition radius, R tr , where the outer 
thin disk transits to the inner ADAF, closes enough to 
the BH. For the convenience of discussion, hereafter we 
refer this type of jets as ergospheric jet from the center 
of a combined disk (EJC). 

Ergospheric jet can also be preferentially launched 
when the BH is surrounded by a slim disk. In this case, 
plasmas enter the BH through a disk-like structure on 
the equatorial plane and most of the region near the hori- 
zon remains nearly vacuum (Figure lb). Therefore, the 
Blandford-Znajek process can efficiently operate. Here- 
after we refer this type of jets as ergospheric jet from the 
center of a slim disk (EJS). The power and the speed of 
the JEC and the EJS are discussed in what follows. 

3.1. Jet Power 

The total extracted power from a rotating BH can be 
computed by P = J R £ r dS, where Rh is the radius of 
the black hole horizon. Similar to Equation (1), we can 
further divide the power into the electromagnetic part, 
P em , and the plasma part, P em , that is, P = P em + 
plasma, where 



dS 



and 



plasma 



Rh 



°plasma ao 



(4) 



(5) 



The mathematical similarity of the idea MHD assump- 
tion and the force-free assumption guarantees that P em is 
identical to the output power of Blandford-Znajek pro- 
cess: 

Pem = ^ 2 F BlRlj 2 C , (6) 

where j = J/Jmax = J/(GM 2 /c), J the angular mo- 
mentum of the BH, £?h the large-scale hole-threading 
field strength, and uj 2 f = ftp(^H — ^f)/^#, the ro- 
tational velocity of the BH. Assuming ujp = 1/2, which 
maximize P em , we have, 



1 1 

324 



2 GJWbH 2 , 2 



(7) 
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Fig. 2. — Thermal equilibria for accretion disks at R/R g = 10 
when Mbh = 1O 9 M0. The solid curves, from left to right, shows 
the ADAF soltuion with a = 0.1, thin disk solution with a = 0.01 
and slim disk solution with a = 0.001. The thin disk solution 
become unstable when the radiation pressure, P ra d, dominates the 
gas pressure, P gas - The filled circle indicates the maximum surface 
density, E, of a thin disk with a fixed a. The maximum X is 
used to evaluate the the maximum large-scale field strength that 
a combined disk can reach (see text for detail). Note that the 
pseudo-Newtonian potential is adopted here. In a full relativistic 
computation, the upper turning point of the "S" curve, which is 
indicated by the filled square an d corresponds to the minimum rh of 
a slim disk, has a value rh > 0.3 (Sadowski 2009; Abramowicz et al. 
120101 ). unlike a large value (rh ^> 1) here. 
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Fig. 3. — Relative large-scale magnetic field strength near the 
BH, £>h, as a function of the accretion-disk type (or, equivalently, 
rh) for a roughly fixed e, according to Equation (3). The increase 
of the field strength with increasing m in a slim disk levels off at 
a saturated value. The resulting magnetically dominated magne- 
tosphere when the BH is surrounded by a combined disk or a slim 
disk, can respectively result in two types of jets, EJC and EJS. See 
text for more details. 

In slim disk case, the accreting plasma, which contributes 
to Ppiasma, stay on the equatorial plane. As a result, 
there is only tiny amount of plasma loading onto the 
large-scale magnetic field lines that thread the hori- 
zon at different latitudes. The jet power of a EJS can 
be therefore estimated by Equation (7), namely that 
^Pejs ~ i^em- On the other hand, because the plasma 
with quasi-spherical geometry near the horizon can be 
loaded onto the field lines and can further reduced the 
power by the amount of P p i aS ma < 0, Equation (7) rep- 
resents the maximum power of a EJC. 

By Equation (3), Bn can be computed if e and E are 
known. Although the determination of e is unclear, it 
should in general has a value 1 > e > 0. Therefore, 
the changing of Bn at different accretion rate should 




10' 10° 
M BH (solar mass) 

Fig. 4. — Maximum power of the EJC as a function of BH masses 
with discrete values of e's. The long-dashed, short-dashed, and 
dash-dotted lines represents the case of £ =0.2, 0.1 and 0.05, re- 
spectively. The power is computed by Equation (5), after the max- 
imum value of £>h is estimated by Equation (3). The spin of the 
BH, j, is assumed to be 0.6. Since the mass loading or a smaller 
spin can further lower the power, the power of EJC will appear in 
the lower part of the se maximum power lin es. The FR 1/ FR II 
division line found by Ledlow & Owen (1996) is represented by the 
solid thick line. Most of FR I (FR II) sources are located above 
(below) this division line. 

mainly be determined by the variation of the surface 
density of the disk, E(ra), because it will be higher vari- 
able than e(m). Here we assume that the large-scale, 
disk-threading fields can be drag ged inward by accre- 
tion flow for all types o f disks (Naravan et al. 2003; 
ISoruit fc Uzdenskvl 120051 : iRothstein fc Lovelacd 12008( 1 
and hence e > is always satisfied. Note that, in a 
slim disk, it is the pressure-driven nature (instead of the 
viscosity) that determines the advection motion of the 
plasma near the BH. Therefore, although the the mag- 
netic Prantdl number is expected to be small due to a 
relative tiny value of a compare to other types of disk 
(Table 1), the field strength near the horizon can be fur- 
ther enhanced when the disk transits form a thin disk to a 
slim disk, by further pushing a large-scale, disk-threading 
field lines inward. 

Figure 2 presents typical BH accretion disc solution 
profiles on the E - rh plane, following t he mo del in 
lAbramowicz et al.l (|1995l ) and IChen et all (jl995l ). The 
solid curves represent the solutions of a ADAF, a thin 
disk and a slim disk with their typical viscosity (see fig- 
ure caption). With the relation of E - rh and Equation 
(3), relative values of Bu as a function of rh are qualita- 
tively presented in Figure 3. The field strength gradually 
increases with increasing rh when the disk is either an 
ADAF or a thin disk because of the increase of E, where 
the radiation-pressure dominated thin disk solution is not 
of interest here because the disk become unstable in that 
case. For a combined disk, Bu is essentially described 
by the surface density of the disk at the transition ra- 
dius, R tr , because the outer thin disk has a much larger 
density than that of the inner ADAF. The decrease of 
R tr with increasing rh results in a rapid growth of B^l 
with rh (see also the Figures la and le of PHC). For a 
slim disk, although Br can be further enhanced due to 
the pushed-in large-scale, disk-threading field lines, the 
increase of Bu is expected to saturate eventually. This is 
because the fields diffuse outward outside the pressure- 
driven region (that is, outside the super-Keplerian part 
of the disk). 
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The transition radius of a combined disk continu- 
ously decreases with increasing rh before the entire disk 
become a thin disk. Previous studies ([Honmal 1996; 
iManmoto et ah 2000) show that R tr can reach down to 
~ 10Rg. For illustration purpose, assuming a modest BH 
spin, j = 0.6, we compute the maximum value of £?h for 
a combined disk by using the maximum E of the (outer) 
thin disk solution at R = 10R g (the filled circle in Fig- 
ure 2) and calculate the resulting maximum power of a 
EJC by equation (3). The jet power for three different 
e's are plotted on the Mbh ~ P plane (Figure 4). Noting 
that that the plasma contribution, P p i a sma < 0, as well 
as a smaller BH spin, will reduce the power of EJC, we 
find these jet powers represent the upper limits. In other 
words, the power of a EJC will appear below these lines. 
On the contrary, for the same parameters (BH mass, BH 
spin, e), sources with EJSs will appear above the lines 
because of a larger £?h (Figure 3) and a negligible P p i aS ma 
in the force- free limit. 

3.2. Jet Speed 

By denoting the quantities of the inflow (u r < 0) and 
outflow (u r > 0) with the subscript "in" and "out", re- 
spectively, the ratio of the jet speeds of a EJC and a EJS 
can be estimated as follows: At large distances, the jet 
Lorentz factor, T, can be defined by 



/i 



(8) 



In addition, near the separation region of the inflow and 
the outflow, the conservation of £ r gives, 



E 



^out^out 



(9) 



Thus, for fixed BH mass and spin, the ratio of the power 
of the EJC and the EJS, %, can be written as, 



X ■ 



T EJC „ (n out < ut ) EJS (B r B^ G 
T EJS ~ (n out < ut ) E JC (B r B^s 



(10) 



by the help of Equation (1), where the superscript 
"EJC" ("EJS") denotes the parameter of EJC (EJS). 

Because both (n ou t< u t) EJS /(^out<ut) EJC < 1 and 
(B r B^ E3C /(B r B^ E3S < 1 hold, we obtain x < 1- 



The characters of EJC and EJS offer the key to an 
understanding of the FR 1/ F R II dichotomy. The FR 
1/ FR II division line found bv iLedlow fc Owenl (|1996[ ). 
which is shown by the thick solid line in Figure 4, is 
roughly consistent with the maximum power that a EJC 
can reach. It is, therefore, reasonable to suppose that the 
division reflects the two types of relativistic jets. That is, 
most FR I (or FR II) galaxies are associated with a EJC 
(or a EJS). Noting that a EJC has a slower speed than 
a EJS, why FR I and FR II galaxies have a edge-darken 
and a edge-brighten morphology, respetively, can be con- 
sistently understood as a result of different jet velocities 
((Bicknell 1985^ 
iGhisellini fc Celottil (120011 ) found that the accretion 



Ghis ellini fc Celottil J2001) 
rate of FR I and FR II galaxies can be separated at 
rh « 0.01. Similar division of the accretion rates is 
also indicated for BL Lacs and radio quasars (which are 
believed to be 'face-on' F R I and FR II g alaxies, re- 
spectively) (jXu et al.ll200i IGhisellini et alJ l2011h . and 
for radio galaxies of 'hig h-excitation' and 'low-excitation' 
(Best & Heckman 2012). The corresponding rh of EJC 
(rh < 0.01) and EJS (rh > 0.01) are consistent with the 
above findings. 

The major drawback in ou r calculation is that, the 
Pseudo- Newtonian potential ([Paczvhskv fc Wiitai ri98Q) 
is adopted when we u se u sed the disk mo del of 
lAbramowicz et al.l (|1995l ) and IChen et al.l (11995ft. We 
should notice a full-relativistic model (e.g. iSadowskll 
(|2009[ )), could significantly change the range of rh of the 
slim disk solution (see captions of Figure 2). Neverthe- 
less, the surface densi ty of the disk solution will not be 
changed dramatically (Sadowski 2009; Abra mowicz et all 



120101 : ISadowski et~atl 1201 ID . We therefore expect that 
the results will be qualitatively unchanged when general- 
relativistic corrections are taken into account. The effect 
of the BH spin on the jet power, which is beyond the 
present model, can be investigated in the future by in- 
corporating general relativistic corrections. 
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007-022-MY3, and the Formosa Program between NSC 
and Consejo Superior de Investigaciones Cientificas in 
Spain administered under the grant NSC100-2923-M- 
007-001-MY3. 
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